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摘要 : 本 文 主要 利用 粗 粒 化 分 子 动力 学 方法 对 不 同 浓度 下 表面 活性 剂 水 溶液 在 couette 流动 过 程 中 胶 束 的 形成 、 速 度 
分 布 与 流 变 性 规律 进行 了 研究 。 结 果 表 明 ， 对 于 相同 浓度 的 表面 活性 剂 溶液 ， 随 着 剪 切 速率 的 增 大 ， 胶 束 完 全 动态 平 
衡 时 间 越 来 越 小 ， 而 随 着 表面 活性 剂 分 子 浓度 的 增加 ， 胶 束 动 态 平 衡 时 间 也 越 来 越 小 ， 形 成 胶 束 越 来 越 快 。 在 剪 切 速 
率 与 溶液 前 切 粘度 关系 方面 ， 保 持 表 面 活性 剂 分 子 浓度 不 变 情 况 下 ， 随 着 剪 切 速率 的 增 大 ， 剪 切 粘 度 逐 渐 减 小 ， 呈 现 
剪 切 稀 化 特征 。 
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A coarse-grain molecular dynamics simulation for surfactant aqueous Solution in 


couette flow 
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Abstract: In this paper, the formation of micelle, velocity distribution, and rheological behavior in 
surfactant aqueous solution under couette flow were studied using coarse-grained molecular dynamics 
simulations. Simulation results show that the dynamic equilibration time becomes shorter with the increase 
of shear rate on the condition of the same surfactant aqueous solution concentration. If the shear rate is a 
constant number, the formation of a micelle was more rapid when the surfactant molecular concentration 
increases in the aqueous solution. When the surfactant molecular concentration keeps constant, the shear 
viscosity becomes smaller with the increase of shear rate, and the shear-thinning characteristic can be 
obtained. 
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0 引言 


随 着 我 们 对 能 源 的 需求 量 越 来 越 大 ， 新 能 源 活性 剂 ， 在 表面 活性 剂 溶液 分 子 动力 学 模拟 方面 
的 开发 以 及 能 源 的 高 效 利用 就 变 得 越 来 越 重 要 ， 的 研究 较 少 ，N.A. Volkov、N. Arai 等 人 的 研究 比 
而 减 阻 作为 一 种 有 效 减 小 管道 内 流动 阻力 的 方 较 有 代表 性 。 其 中 N.A. VolkovI9 对 十 二 烷 基 硫酸 
法 ， 可 以 达到 节能 降 耗 的 目的 。 所 以 减 阻 用 表 钠 dodecyl sulphate, SDS) 溶液 中 胶 束 
活性 剂 的 研究 成 为 目前 的 研究 热点 之 一 44。 而 表 形成 进行 分 子 动力 学 模拟 ， 发 现 NaCl 的 加 入 对 
活性 剂 在 溶液 中 的 微观 结构 是 影响 表面 活性 剂 ww 该 模拟 详细 描述 了 胶 束 
流 变性 的 关键 ， 因 此 深入 了 解 表面 活性 剂 在 溶液 ”的 形状 特性 。N. Arai PUE DPD 方法 研究 了 纳米 
中 的 胶东 形成 和 流 变 性 质 具 有 重要 意义 。 尺度 平行 通道 内 表面 活性 剂 溶液 在 前 切 流 场 中 的 
前 对 于 表面 活性 剂 溶液 的 研究 主要 集中 在 实验 流 变 特性 ， 发 现 通过 对 纳米 平面 通道 内 表面 进行 
测量 阶段 ， 即 通过 宏观 流 变性 的 测量 来 推测 微观 疏 水 - 亲 水 化 学 改 性 , 可 以 使 表面 活性 剂 溶液 实现 
结构 或 者 直接 通过 电镜 对 表面 活性 剂 的 微观 结构 由 剪 切 变 稀 - 剪 切 增 稠 的 改变 。 而 A. 
进行 观察 57。 在 微观 结构 模拟 方面 ,一 般 都 是 针 Sambasivam!?! 则 首次 采用 分 子 动力 学 
对 普通 聚合 物 的 胶 柬 形成 、 断 裂 、 流 变性 等 方面 ( Coarse-grained molecular dynamics, CGMD ) 模 
的 分 子 动力 学 模拟 研究 ， 但 是 具体 到 双亲 性 表 拟 方法 对 CTAC+NaSal 溶液 剪 切 流 场 下 的 棒状 胶 


束 变形 和 断裂 行为 进行 模拟 ， 通 过 能 量 分 析 方 法 (12-60 模型 ， 而 表面 活性 剂 头 基 - 其 它 表 
对 单 体 胶 束 长 度 和 剪 切 流 场 下 胶 束 取向 分 布 进行 。 剂 分 子 尾 基 、 表面 活性 剂 尾 基 与 粗 粒 化 水 分 


了 研究 ， 其 中 得 出 了 CTAC 胶 束 长 度 平均 值 为 间 采 用 LJ9-6 模型 施加 排斥 力 ， A N 
5.4nm 的 结论 ， 该 模拟 的 总 时 间 长 度 为 34s。 在 Hydrophilic 
: e \ 的 定量 研究 字 Head Hydrophobic 
胶 束 形成 时 间 和 胶 束 大 小 的 定量 研究 方面 ， ea Tail 
Wataru Shinoda00 模 拟 了 聚 氧 乙烯 等 几 种 表面 活 i 


性 剂 在 空气 -水 、 油 -水 界面 上 胶 束 形成 过 程 ， 对 
胶 束 大 小 与 胶 束 形成 时 间 进 行 了 定量 计算 。 在 粗 
粒 化 表面 活性 剂 溶液 couette 流动 在 大 剪 切 率 情 一 
况 下 流 变 性 研究 方面 ，J. Castillo-Tejas?b6] -F 7x 一 
烷 基 三 甲 基 氧 化 铵 (CTAC) 表面 活性 剂 溶液 的 

couette 流动 和 poiseuille 流动 进行 了 粗 粒 化 非 平 
衡 分 子 动力 学 模拟 ， 并 且 首 次 对 CTAC 单 体 进行 
了 一 个 亲 水 头 基 和 四 个 疏水 尾 基 的 粗 粒 化 假定 ， 
得 到 了 不 同 表面 活性 剂 浓 度 下 溶液 couette 流动 图 1 表面 活性 剂 与 水 分 子 受 力 示意 

剪 切 速度 与 前 切 应 力 的 关系 ， 剪 切 速 度 与 前 切 粘 Fig.l Surfactant and water CG molecule interaction 

度 的 定量 关系 。 作 为 对 J. Castillo-Tejas 模拟 表面 活性 剂 分 子 头 基 与 相 邻 尾 基 、 尾 基 与 相 
paie E.O. Castrejón- González/?19:3 T 7K 邻 尾 基 之 间 采 用 Harmonic 势能 模型 。 表 面 活性 剂 
-表面 活性 剂 分 子 - 油 溶液 混合 物 在 大 剪 切 率 情 况 溶液 采用 NVT 系 综 ， 固 体 壁面 采用 NVE 系 综 。 
下 couette n n Y 31] 7) 6 BU, 得 到 了 粗 粒 化 表面 活性 剂 溶 液 couette 流动 如 图 2 所 示 。 
随 着 剪 切 速率 的 增加 ， 溶 液 呈 现 由 牛顿 流体 向 剪 
切 变 黎 ， 再 向 剪 切 增 稠 ， 甚 至 最 后 再 向 牛顿 流体 
转变 的 结论 。 
BEME, ERRE REAA, 
究 者 通常 使 用 比较 精细 的 分 子 模型 进行 胶 束 形成 
研究 ， 但 是 这 种 研究 有 两 个 缺点 ， 即 通用 性 差 和 
计算 量 巨大 。 所 以 本 文采 用 粗 粒 化 分 子 动力 学 方 
法 进行 研究 ， 简 化 了 分 子 模拟 所 考虑 的 分 子 极 性 
和 不 同 基 团 之 间 的 受 力 ， 得 到 胶 束 的 形成 时 间 。 
在 对 表面 活性 剂 溶液 的 El, hu d. 


2 


对 表面 活性 剂 的 流 变 性 进行 了 定性 UH. MUT lower fixed wall 

与 目前 宏观 实验 结果 相近 的 结论 ， 加 深 了 对 表 图 2 粗 粒 化 表面 活性 剂 溶液 couette 流动 示意 

活性 剂 溶液 流 变 性 的 认识 。 Fig.2 The surfactant aqueous solution in couette flow 
光 I 1.1 L-J 模型 

1 数值 方法 


L-J (12-6) 势能 表达 式 : 


本 文中 坐标 xyz 方向 的 表面 活性 剂 溶液 模拟 m T2 
域 分 别 为 300 X200 X150 (上 下 壁面 区 域 Y Ur =4e 9 (z) med 
d rV. V. 
坐标 分 别 为 200 -22.50 f0-2.50 -00 )， 其 中 每 * 


个 粗 粒 化 表面 活性 剂 分 子 由 一 个 亲 水 头 基 (直径 。 截断 半径 /2.56 。 
为 4sc ) RU IRSE (直径 为 o ) AR, M 
粒 化 水 分 子 直径 也 为 o 。 其 中 表面 活性 剂 分 子 头 
基 - 表 面 活性 剂 头 基 、 表 面 活性 剂 头 基 - 粗 粒 化 水 
分 子 、 表 面 活性 剂 尾 基 - 表 面 活性 剂 尾 基 、 粗 粒 化 
水 分 子 - 粗 粒 化 水 分 子 之 间 的 作用 均 采 用 LJ 


排斥 势能 函数 CReplusive potential) : 
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T 0.001， 每 组 工 况 运 行 150 732b, £973 5.5ns. 
u= l : 表 1 各 物理 量 无 量 纲 化 公式 
lj Table 1 Dimensionless of parameters 
可 以 用 LJ96 模型 近似 替代 。 3 模拟 结果 与 讨论 


1.2 Harmonic 模型 
Harmonic 势能 函数 表达 式 : 
U,(r,) = 天 (7 -es) 


3.1 胶 束 形成 时 间 分 析 
通过 对 胶 束 形成 过 程 的 可 视 化 观测 《使 用 


物理 量 名 称 无 量 纲 化 公式 
而 0, 通过 Lorentz-Berthelot 法 则 : 
| 长 度 1=1 
c; =0.5(c, +0; ) 计 算得 到 。 Bui 
1.3 NVT 5 NVE 系 综 时 间 ZIXLINC 
NVT 正则 系 综 ( 使 用 Nose-Hoover 人 恒温 条 件 ) 密度 p-pmo^ 
TeT'. 
dr, 剪 切 率 
E =v, pE Èl Em) 
7 * EE 
dv, xS gs f =f e 1 
dt m, ! 压力 p-pe. 
dë e 1 | Di r Ovito 开源 软件 ，Open visualization tool) ， 可 以 
dt Qj|'r 2m, 得 出 couette 流动 剪 切 速率 分 别 为 0.03、0.06、0.1、 


T v, FEAERBISSEE dum, T %15 020 025 六 种 工 况 下 表面 活性 剂 粗 粒 化 


U- o 数 ; 而 Q 定义 为 恒温 质量 , 表达 式 为 D-ZTr，r 为 初始 状态 -小 胺 束 团 状态 -大 胺 束 团 状态 -完全 
为 系统 的 松弛 时 间 。 动态 平衡 状态 。 以 表面 活性 剂 分 子 浓度 0.20, 38 
NVE 系 综 主要 在 壁面 使 用 ， 系 统 的 原子 数 、 切 速率 0.1 为 例 进行 说 明 , 分 别 如 图 3a. 3b, 3c. 
能 量 、 体 积 保持 不 变 。 3d 所 示 。 
1.4 其 它 模拟 细节 


各 个 粗 粒 化 分 子 使 用 Material Studio 进行 初 
始 排 布 ， 排 布 后 的 各 个 粗 粒 化 分 子 坐 标 信息 作为 
ER 模拟 的 初始 条 件 。 此 非 平衡 态 分 子 动力 学 
7 (NEMD) 模拟 中 ， 各 工 况 的 无 量 纲 时 间 步 长 均 
为 
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(a) 初始 状态 0ns (b) 小 胶 束 团 状态 0. 37ns 


(c) 大 胶 束 团 状 态 0. 73ns (d) 完全 动态 平衡 状态 3. 28ns 
图 3 胶 束 动态 平衡 状态 
Np Fig.3 The dynamic equilibrium state of micelles 


N 表面 活性 剂 浓度 分 别 为 0.05、0.15、0.20 下 ”3.2 表面 活性 剂 浓度 与 溶液 速度 分 布 


不 同 剪 切 率 与 胶 束 形成 时 间 的 关系 对 比 ， 如 图 4 研究 表面 活性 剂 分 子 形成 的 胶 束 对 couette 
所 示 。 流动 流 场 的 影响 , 以 couette 流动 前 切 速率 为 0.1、 


胶 束 浓度 分 别 为 0、0.05、0.15、0.20 的 速度 分 布 
ii mcm 为 研究 对 象 进行 说 明 ， 得 到 胶 柬 浓度 与 速度 分 布 
pv 和 各 关系 ， 如 图 5 所 示 。 有 具体 实现 办 法 : 在 Y 方向 上 

( 即 速度 梯度 方向 ) 从 0o 开始 到 200 结束 划分 
出 100 个 等 距离 长 方 体 空间 ， 即 每 个 小 长 方 体 空 
间 的 体积 为 30o X0.20 X156 。 对 每 个 切片 长 
方 体 空间 内 的 粗 粒 化 粒子 (包括 表面 活性 剂 分 子 
和 水 分 子 ) 进行 一 段 时 间 内 的 速度 平均 ， 计 算 系 
ET 


A 


e 
a 


i 


N 
a 


N 


micelle deformation time (ns) 
[^] 


E 
a 


——»——— concentration 0.15 
no surfactant 
———86—— concentration 0.05 


图 4 剪 切 速率 与 胶 束 形成 时 间 ——6—— concentration 0.20 
Fig.4 shear velocity and micelle formation time 
可 以 看 出 ,表面 活性 剂 分 子 浓度 分 别 为 0.05. 
0.15、0.20 情况 下 ， 随 着 couette 流动 上 板 剪 切 速 
率 的 增 大 ( 剪 切 速率 分 别 为 0.03、0.06、0.1、0.15、 
0.20、0.25)， 胶 束 完全 动态 平衡 时 间 越 来 越 小 。 
而 在 相同 剪 切 速率 情况 下 , 随 着 胶 束 浓 度 的 增加 ， 
胶 束 动态 平衡 时 间 也 越 来 越 小 ， 即 胶 束 形成 速度 T 
也 越 来 越 快 。 Y coordinate 


0.1 0.15 2 
shear rate (reduced units) 
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图 5 胶东 浓度 与 空间 平均 速度 分 布 
Fig.5 The relationship between micelle concentration and 
the spatial mean velocity distribution 
从 图 5 可 以 看 出 , TEBEUI SERES 0.1 情况 下 ， 
水 溶剂 《表面 活性 剂 浓度 为 0， 图 示 中 黑 线 表示 ) 
的 速度 分 布 基本 是 单调 递增 的 ， 成 斜率 为 1 BUE 
性 分 布 。 在 加 入 粗 粒 化 表面 活性 剂 分 子 以 后 ， 在 
Y Jj i] 0-11 0 区 域 平均 速度 大 于 水 溶剂 ， 而 
11-20 c: 区 域 速度 则 小 于 纯 水 。 至 于 是 否 表 面 活 性 
剂 浓 度 越 小 ， 速 度 分 布 偏离 纯 水 越 小 ， 模 拟 中 的 
数据 不 予 支持 ， 因 为 表面 活性 剂 分 子 浓度 为 0.15 
时 在 11-20 0 区 域 与 浓度 为 0.20 时 相 比 更 加 偏离 
纯 水 的 速度 分 布 。 所 以 可 以 得 出 一 个 结论 ， 表 面 
活性 剂 浓 度 大 小 与 改变 couette 流动 速度 分 布 的 


幅度 大 小 没有 必然 关系 。 
33 剪 切 速率 与 剪 切 粘度 关系 
剪 切 粘度 由 前 切 应 力 与 前 切 速率 得 到 ， 计 算 
公式 为 ; 
- Py 
TK 
其 中 为 剪 切 粘度 ， 已 ,为 剪 切 应 力 ， 双 为 前 


切 速率 ， EA A T Deme 

无 量 纲 温度 T=1. 35. 76 EC ApE 9E O7] RIIA 
0.01、0.025、0.05、0.1、0.15、0.2、0.25 I], Bj 
切 速率 与 剪 切 应 力 、 剪 切 粘 度 的 关系 ， 分 别 如 图 
6、 图 7 所 示 。 


——8——- concentration 0.05 
——-8——— concentration 0.10 
——8——— concentration 0.15 
——*e— concentration 0.20 


0 0.05 0.1 0.15 02 0.25 
shear rate 
图 6 剪 切 速率 与 剪 切 应 力 关 系 


Fig.6 The relationship between shear rate and shear stress 


———8——- concentration 0.05 
——-8——— concentration 0.10 
——8——— concentration 0.15 
-一 全 一 concentration 0.20 


一 一 


shear viscosity 
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图 7 剪 切 速率 与 前 切 粘 度 关 系 
Fig.7 The relationship between Shear rate and Shear 
viscosity 

从 图 6 可 以 得 出 如 下 规律 随 着 剪 切 速率 的 
增 大 ， 剪 切 应 力 随 之 增加 。 而 随 着 粗 粒 化 表面 活 
性 剂 分 子 浓 度 增 大 ， 剪 切 应 力也 越 大 。 

从 图 7 可 以 看 出 ， 除 了 剪 切 速率 为 0.01 和 
0. 025 If DI WREE AP f ERAL, 没有 呈现 出 剪 
切 稀 化 的 特征 之 外 ， 之 后 随 着 剪 切 速率 的 增 大 ， 
剪 切 粘度 逐渐 减 小 ， 呈 现 剪 切 稀 化 特征 。 剪 切 速 
率 太 小 (=0.03) ， 粗 粒 化 分 子 的 热 运 动 干扰 过 
大 ， 导 致 系统 内 各 层 粗 粒 化 分 子 不 能 完全 做 有 序 
的 couette 流动 ,所 以 只 取 0. 05-0. 25 有 效 剪 切 速 
率 区 域 进行 分 析 。 同 时 可 以 看 出 粗 粒 化 表面 活性 
剂 分 子 浓度 与 剪 切 粘 度 之 间 呈 现 正 相 关 ， 即 表面 
活性 剂 浓 度 越 大 ， 剪 切 粘 度 越 大 。 


— 


4 结论 w 


本 文 主要 对 不 同 浓度 的 表面 活性 剂 溶液 
couette 流动 进行 了 非 平衡 分 子 动力 学 模拟 , 得 出 
couette 流动 中 剪 切 速 率 与 粗 粒 化 表面 活性 剂 分 
子 胶 束 形成 时 间 、 表 面 活性 剂 浓度 与 游 液 速度 分 
布 、 剪 切 速率 与 剪 切 应 力 、 剪 切 速率 与 前 切 粘度 
的 关系 。 发 现 表面 活性 剂 浓 度 和 剪 切 速率 与 胶 柬 
形成 时 间 负 相关 ， 表 面 活 性 剂 浓 度 与 溶液 速度 分 


布 无 必然 联系 。 在 浓度 一 定 情况 下 ， 剪 切 速率 越 
大 ， 剪 切 应 力 越 大 ， 剪 切 粘 度 越 小 ， 呈 现 剪 切 稀 


化 特征 。 通 过 对 表面 活性 剂 溶液 中 这 些 特定 参数 
的 分 析 研 究 ， 可 以 更 深入 地 了 解 表面 活性 剂 溶液 
胶 束 的 形成 以 及 胶 束 的 流 变 特性 。 
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